1 



Effective Potential for Polyakov Loops in Lattice QCD 

Y. Nemoto a [RBC collaboration] * 

a RIKEN BNL Research Center, Brookhaven National Laboratory, NY 11973, USA 



(N 
O 

o 

(N 
•*-> 
O 

O 



> 

00 

o 
o 

o 



X 



Toward the derivation of an effective theory for Polyakov loops in lattice QCD, we examine Polyakov loop 
correlation functions using the multi-level algorithm which was recently developed by Luscher and Weisz. 



1. Introduction 

Precise measurement of Polyakov loop correla- 
tion is important for finite temperature physics. 
First of all heavy quark potentials at finite 
temperature are obtained by measuring two- 
Polyakov-loop correlation functions. It is already 
well-known that we can obtain the QQ potential 
below and above the critical temperature (Tc) 
from this measurement. Due to the nonpertur- 
bative effects, however, there are still less-known 
quantities at intermediate temperature. For ex- 
ample the QQ potential at finite temperature is 
generally written as V/T = -e/(RT) d exp(-(j,R), 
where T is temperature, R the distance between 
quarks, and e, d, n are parameters. At very high 
temperature, perturbation theory predicts d = 2, 
while just above Tc, lattice simulation gives the 
deviation from d = 2 fjj . In relation to such a non- 
perturbative behavior at finite temperature, our 
main interest here is to derive an effective theory 
for Polyakov loops from lattice QCD. Recently 
new Polyakov loop model is proposed^) regard- 
ing the deconfining phase transition in SU(N) 
gauge theory. It conjectures a relationship be- 
tween the Polyakov loops and the pressure. The 
usual Polyakov loop l\ transforms under a global 
Z(N) transformation as l\ — > e ll ^l\, with <fi = 
2irj/N(j = 1...N — 1). This means that l x trans- 
forms as a field with charge one. On the other 
hand one can consider Polyakov loops with higher 
charge under a global Z(N) transformation. For 
example the charge two Polyakov loop is defined 
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where L is an ordered link product in the tem- 
poral direction, L(x) = Pexp(ig Aq(x, r)dr). 
The Polyakov loops with higher charge might af- 
fect non-universal behavior in the model. There- 
fore it is interesting to measure relevant kind 
of new Polyakov loop correlation on the lattice 
without any phenomenological assumptions. In 
this case, however, we need the Polyakov loops 
with relatively large length in the temporal di- 
rection, which might be inaccessible in the ordi- 
nary computation of the Polyakov loops. Thus 
we are led to try the new algorithm recently de- 
veloped by Luscher and Weisz called the multi- 
level method to measure various Polyakov loop 
correlation functions. This method can measure 
quite precisely for the Polyakov loop correlation 
with the large temporal length and/or large dis- 
tance between the Polyakov loops as shown be- 
low. 

2. Precise measurement of Polyakov loop 
correlation 

One way to get the Polyakov loop correlation 
functions as precisely as possible is the method of 
improved estimator. When we measure Polyakov 
loops, instead of a time-like link variable itself, 
by using some average of the time-like links, we 
can reduce the error of the Polyakov loop without 
affecting the average. The multi-level algorithm, 
which is used here, is more efficient in measuring 
the Polyakov loop correlation. We explain this 
algorithm briefly in the following. Let us consider 
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In this method there are essentially two new 
parameters, which we call n ms and n up . n ms is 
the number of measurement to be made for the 
time slice average 



[T(t ,r)T(t + a,r)] 

n ' * 



(6) 



Figure 1. A two-time-like link operator (denoted 
by thick lines, left figure) and two-Polyakov-loop 
correlation with the multi- level algorithm, eqs.(Q) 
or (||) (right). 



the two-Polyakov-loop correlation, (P(n)P^ (n + 
r)}, where n is a spatial coordinate and r is the 
distance between two quarks. First we define a 
two-time-likc link operator, 



T(t a ,r) afil s = U4(to,n) a /3Ul(t ,n + r) l5 . 



(2) 



The essence of the multi-level algorithm is similar 
to the multi-hit one, i.e., is to take an average of 
T instead of T itself in the sublattice of the time 
slice, 
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D{U] s T(t ,r)e~ s M°. (3) 



This is schematically shown in Fig.l. Here s 
denotes the sublattice which includes the oper- 
ator T and is shown as the shaded area in Fig.l. 
The Polyakov loop correlation function is then 
expressed as 

{P(n)P\n + r)) 
= ([T(0,r)T(a,r)][T(2a,r)T(3a,r)] 

■■■[T(T-2a ir )T(T-a,r)]) (4) 

for single layer or 

(P(n)P\n + r)) 
= ([[T(0,r)T(a,r)][T(2a,r)T(3a,r)]] 

••■[••• [T(T - 2a, r)T(T — a, r)]]) (5) 

for two layers. These equations are also schemat- 
ically shown in Fig.l. 



and n up is the number of time slice updates be- 
tween measurement. We determine n ms so that 
the s/n ratio is equal to unity. Then the statistical 
errors are reduced exponentially as the distance 
between the Polyakov loops increases |3j. 

We show some test results of the algorithm in 
the following. All the simulations are done with 
the lattice size 12 4 at (3 — 5.7. We use the stan- 
dard Wilson action and no gauge fixing. Fig. 2 
shows the n ms dependence of a two-Polyakov-loop 
correlation function at a certain gauge configu- 
ration. In the figure, the spatial distance be- 
tween the two Polyakov loops is set to 4 in the 
lattice unit and n up = 1 is fixed. In the single 
layer case, the relative computation time in the 
horizontal line is proportional to n ms , i.e., the 
leftmost point corresponds to n ms ~ 10 and the 
rightmost n ms = 1000. In the two layer case, 
we have chosen the parameters n ms = 10 and 
n up = 16 for the outer layer and n ms = 5 ~ 150 
and n up = 1 for the inner layer. Usual measure- 
ment of the Polyakov loop correlation without 
any improvement corresponds to n ms = 1. The 
graph shows that the values approach the precise 
one as n ms increases. We have confirmed that 
our calculation agrees with that of Q within the 
range of the statistical error. 

We have also calculated three-Polyakov-loop 
correlation shown in Fig. 3. As an example of the 
result we show the n ms dependence of the quan- 
tity {\P(n + 4x)P(n + 4y)P(n + 4z)\) in Fig.4. We 
see the validity of the algorithm even in this case 
by increasing the value of 

Finally we comment on our trial computation 
of a variant of the multi-level method. Measure- 
ment of many-point Polyakov loop correlation in 
the above method spends much memory. One 
possible solution is to apply a method which is 
used in the calculation of quark propagators of 
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Figure 2. Two-Polyakov-loop correlation with the 
distance r = 4a. 



disconnected diagrams. Let us consider for sim- 
plicity the two-Polyakov-loop correlation. The 
modified method is to take summation over space 
of the two-time-link operators before taking the 
tensor product of them, 



(P(n)P\n + r)) 
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■{[[K(0)][K(2a)]\ 



N con { 

■■■[[K{T-Aa)][K{T-2a)]]}, 
K(to) = T(t ,r)T(t + a,r). 



(7) 
(8) 
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If the number of configuration and n ms are suffi- 
ciently large, this method gives the same results 
as the conventional multi-level method, because 
the gauge non-invariant terms cancel out. Our 
trial simulation is, however, much noisier than the 
conventional one due to the finite number of con- 
figuration and n ms . So unfortunately our present 
results imply that this method is less effective. 

In summary, toward the construction of the ef- 
fective theory for Polyakov loops, we have tried 
a new algorithm called the multi- level method ||. 
We have reproduced the two-Polyakov-loop cor- 
relation by H and confirmed the validity for the 
computation of the three-Polyakov-loop correla- 
tion. 



Figure 3. Schematic figure of the three-Polyakov- 
loop correlation. 



rx=ry=rz=4 



10" 



N 10"' 
+ 



+ 



10 



£L 10"' 
x 

+ 



10" 



• 1 layer 



1000 2000 

nms 



3000 



Figure 4. n ms dependence of the three-Polyakov- 
loop correlation function. 
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